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Abstract 

The ground state of a system of N impenetrable hard core quantum particles in a 
1-D box is analyzed by using a new scheme applied recently to study a similar system 
of two such particles [Centl. Eur. J. Phys., 2(4), 709 (2004)]- Accordingly, each 
particle of the system behaves like an independent entity represented by a macro- 
orbital, -a kind of pair waveform identical to that of a pair of particles moving with 
{q, —q) momenta at their center of mass which may have any momentum K in the 
laboratory frame. It concludes: (i) < A6{x) >= 0, (ii) < x >> A/2 and (iii) q > 
Qoi= T^/d) (with d = L/N being the average nearest neighbor distance), etc. While 
all bosons in their ground state have q = qo and K = fermions have q = qo with 
different K ranging between and K = Kp (the Fermi wave vector). Independent 
of their bosonic or fermionic nature, all particles in the ground state define a close 
packed arrangement of their equal size wave packets representing an ordered state 
in phase ((/>— )space with A0 = 2mi (with n = 1,2,3, ...), < x >= A/2 = d, and 
q = qo- As such our approach uses greatly simplified mathematical formulation and 
renders a visibly clear picture of the low energy states of the systems and its results 
supplement earlier studies in providing their complete understanding. 
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1. Introduction 

Wave mechanics of a 1-D system of identical particles interacting through hard core 
(HC) potential of zero or non-zero range has been a subject of great importance for the 
last seven decades [1-21]. The development of the subject has been facilitated by the fact 
(demonstrated by first few studies by Tonks [1], Girardeau [2,3], and Lieb and Liniger 
[4]) that the problem is exactly solvable. Identical conclusion has been made later [22, 
23] for a similar system of particles interacting through inverse square potential. These 
developments have been elegantly reviewed by Popov [24], Korepin, et al [25], Mattis [26] 
and most recently by Sutherland [27]. As evident from the number of such studies [28-44], 
interest in the subject is renewed recently after the use of Bose Einstein condensation 
(BEG) in dilute gases in experimental realization of many body systems [45-54] that can 
be described by 1-D Hamiltonian with 5— potential interaction [55, 56]. The scheme of 
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theoretical analysis of these systems uses Bethe ansatz [57] for N body wave function 
with bosonic/fermionic symmetry (as the case demands) and suitable (periodic/cyclic) 
boundary conditions. Proposed initially by Lieb and Liniger [4], properties of such a 
system can be determined theoretically by introducing a dimcnsionless parameter, 7 = 
'mgiD/nh^ (with m,-the mass of a particle, n, -the number density of particles and Qid, -the 
1-D coupling constant related to 3-D s— wave scattering length a). One identifies Tonks- 
Girardeau (TG) regime of strong coupling (large Qid) and low density (n) by 7 ~ cxd, 
Bogohuobov regime [56] (or decoherent regime [53]) of nearly free bosons by 7 0, and 
Lieb and Liniger (LL) regime of weak coupling (low qid) and high density by intermediate 
values {viz., 00 > 7 > 0). One also defines Gross-Pitaevskii (GP) regime [53] by 1 > 7 > 0. 
Apparently, the behavior of a system of given giD, can be studied in its different regimes 
simply by changing n from to 00. One of the important conclusions of these studies 
reveals that the physics of a bosonic system characterized by 7 f=:i 00 resembles with that 
of free fermions. 

In this paper, we analyze the ground state (G-state) of a 1-D system of N HC particles 
by using a shghtly different scheme. The excitations of the system would be investigated 
in our forthcoming article. Our scheme uses the wave mechanics of a pair of HC particles 
{identified as the basic unit of the system) as its basis; it has been investigated in our recent 
paper [58] for two HC particles trapped in a 1-D box [58]. Our approach renders exact 
solutions. In view of the equivalence of impenetrable 5— repulsion with HC interaction of 
non-zero range [i.e., Vhc{x), defined by Vnci^ < cr) — 00 and Vhc{x > cr) = with a 
being the HC diameter of a particle] demonstrated by Huang [59] and physically argued 
in [58], one may find that our results for 5-repulsion can be applied to particles of any a, 
particularly, in the low energy states of A/2 > a. 

Our results differ from those of extensively used scheme(s) of others [24-27], particu- 
larly, for the low energy states where the well known consequences of wave particle duality, 
viz., inter-particle phase correlation, zero-point repulsion, quantum size, etc. (defined and 
discussed in Section-3 for their better understanding), become important. Our scheme 
provides a mathematically simple, unified, and clear picture of 1-D systems. As such 
one may find that the present study supplements earlier studies in providing a better 
understanding of our 1-D systems and its basic results resemble with those of : (i) a 
single particle in 1-D box [60, 61], two HC particles in a similar box ([58], and (iii) long 
awaited microscopic theory of a 3-D system of interacting bosons [62-64] which explains 
the properties of liquid ^ife (a well known representative of such systems) with unmatched 
accuracy, simplicity and clarity. This demonstrates the accuracy and usefulness of our 
scheme. As reported briefiy in [65, 66], our scheme provides a theoretical framework 
that unifies the physics of widely different systems (including low dimensional systems) 
of interacting bosons and fermions, atomic nucleus, newly discovered BEG state, etc. and 
helps in revealing the basic foundations of the microscopic theory of superconductivity. 

The paper has been arranged as follows. While the important aspects of our N- 
body system are described in Section-2, the important inferences of our study of two HC 
particles in 1-D box, which serves the basis of this work, are summed up in Section-3. 
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This is followed by a detailed analysis of the ground state of the system in Section-4 and 
concluding remarks in Section- 5. 



2. Important Aspects of N-Body System 

The hamiltonian of a 1-D system of N particles interacting through a two body im- 
penetrable 5-potential can be expressed as 

where A, representing the strength of Dirac delta function {S{xij)), is such that A — > oo 
with Xij — > 0, while other notations have usual meaning. Possible interactions such as 
spin-spin interaction are not included in Eqn. 1 by presuming that these can be treated 
perturbatively. 

Analysing the basic nature of the possible dynamics of two particles (say PI and P2) 
in the system, we note that: (i) PI and P2 encounter 5-potcntial only when they suffer 
a collision, (ii) in a two body collision, they simply exchange their momenta, ki and 
k2, while in a many body collision, where they also collide simultaneously with other 
particle(s), they could be identified to jump from their state of ki and k2 (or their relative 
momentum A; = k2 — k2 and center of mass (CM) momentum K = k2 + ^2) to that of new 
momenta k'l and k'2 (or k' and K'), and (iii) between two collisions each of them has free 
particle motion. Evidently, a state of two particles even in the presence of other particles 
of our system can be identified characteristically with the state of only two HC particles 
in 1-D space. In other words a pair of particles forms the basic unit of our system and 
this agrees with the fact that particles in the system interact through a pair potential. 

3. Dynamics of two HC particles (PI and P2) 

We note that the dynamics of PI and P2 can be described by the Schrodinger equation 
for H{2) [Eqn. 1 with = 2] expressed in the CM coordinate system as 

with 

^(x, X) = i^kix) exp [iiKX)] (3) 

where ipkix) and exp [i{KX)], respectively, define the relative and CM motions of PI and 
P2. We have 

X = X2 — xi and A; = k2 — ki = 2q, (4) 

X = (xi+X2)/2 and X = ki + k2 (5) 

with X and k, respectively, describing the relative position and relative momentum of PI 
and P2, while X and K, similarly, defining the position and momentum of their CM. 
Without loss of generality, one may also define 

ki^-q+— and k2 = g y . (6) 
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It is evident that ipk{x) is a solution of 



-—-g^ + A6{x)\Mx) = EkMx) (7) 



X 



where = E — K"^ / Am. With a view to find ipk{x), we note that PI and P2 experience 
zero interaction at x 7^ where each of them can be described by separate plane waves such 
as Mki(xj,t) = exp (ikjXj) exp (— (assumed to have unit normalization), however, 
for the possible superposition of these waves, their quantum state is better expressed by 



*(xi,X2,t)=^ = l/V2.[Mki(xi,t)Mk2(x2,t) ±Mk2(xi,t)Mki(x2,t)]. (8) 

Since the state function of two impenetrable HC particles must vanish at Xi = X2, 
^(xi, X2, t)"*" (having +ve symmetry for their exchange) does not represent the desired 
function, while ^(xi,X2,t)~ of —ve symmetry has no such problem. We addressed this 
problem in our recent study [58] of the wave mechanics of two impenetrable HC particles 
in 1-D box and used another method to obtain the right wave function of +ve symmetry. 
Evidently, using our results of [58], we express a state of two HC fermions/bosons (in the 
CM coordinates) by 

C(x, X, t)^ = Ck{x, t)^ exp [i{KX)] exp [-i{Ek + EK)t/h] (9) 

with 

(^[x, t)~ = V2 sin {kx/2) exp [-iEkt/Ti] (10) 
of fermionic symmetry and 

Ck{x, t)+ = V2sm {\kx\/2) exp [-iEkt/H] (11) 

of bosonic symmetry. Note that ({x,X,t)~ is a form of ^(xi,X2,t)~ (Eqn. 8) expressed 
in CM coordinates. Although, two impenetrable HC particles in 1-D are expected to 
retain the order of their locations {e.g. < Xi < X2 < L for two particles in a 1-D box 
with infinite potential walls located at x = and x = L) but Eqns. 9-11 remain valid for 
describing the state of PI and P2 since the two particles during their collision do exchange 
ki and k2 which is equivalent to exchanging sides because the particles are identical. In 
fact in a state of wave mechanical superposition of two identical particles, one has no 
means to ascertain whether the particles exchanged their momenta or their positions and 
we may, justifiably, analyze Eqns. 9-11 to find important aspects of the wave mechanics 
two HC particles as follows. 

(i) SMW state : ({x, X, t)"^ represents a kind of standing matter wave (SMW) which 
modulates the probability, \({x, X,t)'^\'^ = |Cfc(.x, t)''=p, of finding the two HC bosons (or 
fermions) at their relative phase position = kx in 0— space. However, the equality, 
|CA;(a;, t)~p = \(k{x,t)~^\'^, renders an important fact that the relative configuration of 
PI and P2 in (^^(a;, t)^ states is independent of their fermionic or bosonic symmetry. 
Although, dependent terms in Eqns. 9-11 have their own importance and they help in 
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having a better understanding of the SMW nature of ^(x, X, t)^, nevertheless these can 
be dropped for their no impact on our impending analysis. 

(ii) Symmetry of relative motion : The relative motion of PI and P2 maintains a center 
of symmetry at their CM and this implies that 

xcm(1) = -xcm(2) = 2;/2 and kcM(l) = -kcM(2) = g (12) 

where xcm and kcM, respectively, represent the position and momentum of a particle with 
reference to the CM of PI and P2. Evidently, two particles in C,k{x,t)^ state have equal 
and opposite momenta (g, —q) with respect to their CM whose motion in the laboratory 
frame can be defined by a plane wave {ex^{iKX)) of momentum K and this agrees with 
what we learn from Eqn. 6. 

(iii) MS and SS states : We note that ({x, X, t)^ resulting from the superposition of two 
plane waves of ki and k2 could be identified to represent the mutual superposition (MS) 
of PI and P2 because ({x, X, t)"^ is, basically, an eigenstate of the energy operators of the 
relative and CM motions of two particles rather than of individual particle. However, it 
could also be identified as a state of the self superposition (SS) of individual particle (PI 
or P2) by identifying C(x, X, t)='= as a macro-orbital {cf. Point- vii below) on the basis of 
the following analysis. To understand the meaning of self superposition one may track 
down the motions of PI and P2 separately and find that each of them (say PI), after 
its collision with P2, has superposition of its pre-collision plane wave ■Ukj(xi,t) with its 
post-collision plane wave Mk'^(xi,t) (with its new momentum k'^ = k2 because during 
the collision PI exchanges its momentum with P2 and vice versa); similar observation 
applies to P2. As such (kixjt)"^ part of our SS state is not different from the state of a 
particle trapped in 1-D box. However, since PI and P2 arc identical particles we have 
no means to identify whether the two particles have their mutual superposition or a self 
superposition of individual particle. Consequently, C{x,X,t)^ could be used identically 
to represent either the MS state of PI and P2 or the SS state of either particle. It may be 
mentioned that two plane waves of two particles can, in principle, have their superposition 
independent of their relative position x and relative momentum k, however, as shown by 
experiments the wave nature of particles seem to influence the behavior of a system of 
particles only when their A = 27r/g compares with their x which implies that an effective 
superposition of two particles leading to their MS/SS state exists only in low energy state 
with A ~ X. 

(iv) Characteristics of relative motion : Cfc(x,t)^ has a series of antinodes of size A/2 
between different nodes x — ±nA/2 (n = 0,1,2,3,...) where A = 27r/g represents 
the de Broglie wave length related to the relative motion of PI and P2. Evidently, two 
particles can be confined to the shortest possible space of A size without disturbing their 
Cfe(a;, t)^ state (say, by placing two impenetrable potential walls at the nodal points, 
X — ±A/2 and for such a confinement we have [58] < x >= I' / 1 = A/2 with integrals 
/' =< C,k{x,t)^\x\C,k{x,t)'^ > and / C,k{x,t)^\C,k{x,t)^ > performed between a; = 
(the minimum possible x) to x — X (the maximum possible x) . Evaluating similar results 
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for < > (as shown in [58]), we find that a state two HC particles in free space is 
characterized by 

<x>> A/2, and <</>>> 2tt. (13) 

This imphes that any experiment, that keeps track of the relative dynamics of two HC 
particles, would find that PI and P2 never reach closer than < x >„— A/2 (shortest 
possible < X >) and < >o= Stt (shortest possible < >). Similarly, Eqns. 12 and 
13 reveal that the position of PI and P2 as seen from their CM can not be closer than 
< xcm(1) >o= —A/4 and < xcm(2) >o= A/4 or vice versa. Finally, we also find that 

<A6{x)>^\Ck{x)^\l=o-0, (14) 

which has been analyzed for its general validity in [58] (see Appendix-A of arXiv.org/ quant- 
ph/0603233) which clearly shows that < AS{x) >— is vahd for all physically relevant 
situations rendering 

< H{2) >= {h^/Am){K^ + e) = {h'/2m){kl + k^). (15) 

(v) Quantum size of a particle : In view of Eqn. 13, it is evident that two particles by 
themselve do not assume a state of < a; >< A/2 (with A = 27r/g) which indicates that 
each of them exclusively occupies A/2 space; being identical particles with equal they 
are expected to share A space equally. Hence as used in this paper we identify A/2 as 
the effective size of a particle particularly for low energy states oi q < njo and name 
it as quantum size. For a better understanding of this meaning, one may consider the 
wave associated with PI as a probe to scale the size of P2 or vice versa and apply the 
principle of image resolution to find that PI (or P2) would be able to resolve the a size 
of P2 (or PI) only if A/2 < a {i.e. q > n/a) implying that in such a case they see each 
other as particles of size a. However, in case A/2 > a, they can not resolve a and would 
see each other as the objects of size A/2 limited by their capacity to resolve. As such 
the effective size of PI and P2 is a g— dependent quantity for X/2 > a {or q < n/a) and 

independent quantity for A/2 < a (or q > n/a). Evidently, the use of quantum size 
for A/2 is justified to distinguish the two situations. This meaning seems to qualitatively 
agree with the meaning of '^quantum spread of a particle^' as used by Huang [59]. However, 
on quantitative scale it seems to represent the minimum possible quantum spread {limited 
by the uncertainty principle) of a particle of momentum q but it, evidently, differs from 
the meaning of "quantum size" as used in relation to the quantum size effects on the 
properties of thin films and small clusters of atoms [67]. The fact that a particle in its 
ground state in a box has a spread of A/2 = d {the size of the box) also indicates that 
the particle exclusively occupies A/2 space because any attempt to decrease this space 
size (say by reducing d) would push the particle to have new momentum/energy but once 
again the space occupied would be a new A/2. 

(vi) Zero point repulsion : Since two HC particles do not assume a state of<a;><A/2 
(c/., Eqn. 13), one would obviously like to identify the responsible force. To this effect 
we note that < x >< A/2 represents an overlap two particles of effective size A/2 and 
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this imphes that particles in this state have higher energy in comparison to the state 
of < X >> A/2. Evidently, the two particles in the state of < x >< A/2 experience 
a kind of mutual repulsion (or the zero point repulsion) unless they assumes a state of 
< X >> A/2. Thus A/2 represents the effective size of a HC particle oi q < n/a and the 
range of zero-point repulsion. 

(vii) Macro-orbital representation of a particle : Although, two particles in C(x, X)^ 
states can be identified to have inter-particle phase correlation ((?(^) = |C(a;, = 
|sin(0/2)p) which can keep them locked at < >= 2n7i {n = 1,2,...) in 0— space, 
nevertheless they do not form a bound pair in a;— space because they either experience 
mutual repulsion when <a;><A/2orno force when < x >> A/2. Evidently, each 
of them in ({xjX)"^ state represents an independent entity described by a separate pair 
waveform, say, ^(a:(j), X(j)) = ({x,X); the subscript (i) refers to i—th particle. To dis- 
tinguish ^(,T(j), X(j)) from ({x,X), we propose to call the former a macro-orbital [58] and 
obtain the same by replacing x, X, k and K in ({x,X)) by and /^(j), 
respectively. This renders 

C(x(i), X(i)) = BCq^.^{x^i)) exp [i(X(i)X(i))] (15) 

where B is the normalization constant and Cg(i)(^(i)) is that part of a macro-orbital which 
does not overlap with similar part of other macro-orbital. In view of what has been argued 
in point (iii) above, the macro-orbital representation is consistent with MS/SS state of 
PI and P2. 

(viii) Macro-orbital as an eigenfuntion : A macro-orbital (^(a;(j), X(j))) is a derived form 
of a wavefunction which identifies a particle to have two different motions (c/. Eqn. 15) 
: the g-motion of energy E{qi^i)) = q'^--^ / 2m = ^^/c^^-)/8m which decides its quantum size 
A(j)/2 = 7i/q(i) and i^-motion of energy E[K{i)) = h'^K'^-^/Sm which represents a kind 
of free motion of the particle. Evidently, this form does not fit, as a solution, with the 
form of Schrodingcr equation (Eqn. 2). However, the corresponding hamiltonian can be 
rearranged to obtain Eqn. 2 in a form with which ^(a;(j), X(j)) is compatible and to this 
end we define 

hi--^^, M^) = ^^^^^ and H{2)^h{l) + hi2) (16) 

with i = 1 or 2 for a system of iV = 2, /?.Ar+i = hi, hi being the kinetic energy operator of 
i-th particle in unpaired format of PI and P2, and h{i) is the same in their paired format. 
We have 

Smdl,,. 2m ^ ^ 

which is such that 



hm^i^),Xi^)) = mqh) + E{Kio)/mHi)^Xit))- (is) 



4. N-Particle System and Macro-orbital Representation 
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(i). Rearrangement of H (N) : For the consistency of macro-orbital representation, we use 
Eqn. 16 to rearrange 



N 



H{N) = ^hi + Y.AS{x,,] 



i>j 



as 



or 



N 



H{N) = Y.h{t) + j:A6{ 



N 



i>j 



8m dx 



i>j 



(19) 
(20) 
(21) 



(ii). State Functions : Using N macro-orbitals for N particles, we construct a state 
function of the system by following standard procedure. Assuming that ^„ represents 
n-th state of energy E^, we have 



with 



and 



0n(g)-n,t,(2/L)^/2sin[g(,)X(,)], 



(22) 



(23) 



(24) 



Since each particle occupies a cage between its neighbors, in Eqn. 23 should be 
integer multiple of 7r/d(j) (d(j) = size of the cage); to a good approximation, q^^^ can also 
be assumed to be an integer multiple of tt/L if L is macroscopically large. In Eqn. 24, we 

— 1 /2 

have B = {N\L^) as the normalization constant, i^'(i) = Nvr/L (with N = 1, 2, 3, ...) 
and X]p(=tl)^ as the sum of A^! products of N plane waves obtained by permuting different 
Ar(j) over N particles with P being the number of performed permutations. While the use 
of (+1)^ or (—1)^ depends on the bosonic or fermionic nature of particles, the condition 
> T^/d{i) follows from the HC nature of particles (bosons and fermions alike). Since 
different particles can, in principle, have different their permutation over particles 
renders A^! different 0n(Q') and, therefore, similar number of Consequently, the general 
form of a state function that should reveal the physics of the system should be expressed 

by 



(AT!) 



TV! 
7 = 1 



(25) 



(iii). Energy Eigenvalues : The energy eigenvalue En can be expressed by 





Ho{N) + Ei<j A6{x,j)\<!>n > _ < $n 


Ef M^)l<^>n> 






$ > 



(26) 
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which renders 



N 



En=i: 



8m 2m 



(27) 



since < Ab{xij) > vanishes (c/., Eqn. 14) for every pair of HC particles. One can also 
conclude < AS{xij) >= o by using A*" body wave function and the detailed derivations, 
as used to conclude similar result for a 3-D N body system [63] . 

(iv). G-state Energy : While has zero value for all bosons in their G-state, its different 
values ranging between and Kp (Fermi wave vector) in the G-state of a fermionic system 
follow Fermi-distribution. To fix possible values of ^(i) for which Eg is minimum, we 
note that each particle satisfies < d{i) which implies that A(j)/2 space exclusively 

belongs to a HC particle of momentum ; this agrees with the volume excluded condition 
envisaged by Kleban [68] for HC particles in 3-D systems. Since each particle in the G- 
state has lowest energy or largest possible A(j)/2, the net energy of g-motion can be 
expressed as 

= X] i — 12~ ^^^^ ^(i) ^ (constant) (28) 

by assuming that particles occupy cages of different size d^iy Simple algebra reveals that 
Eo has its minimum value for ci(i) = d(2) — ci(3)... = c?(Ar) = d. Obviously, for a bosonic 
system we have 

which does not differ from similar result for 3-D bosonic system [64]. Accounting for an 
additional contribution of X-motions to the G-state energy of a fermionic system, we have 

Ef/'"'^ = Ne,+ ^eK, (30) 

where is the total energy of X-motions at T = in a system of non-interacting fermions. 
The factor (|) accounts for the fact that each particle of our system in its macro-orbital 
representation serves as an entity of mass 4m for its K motion. For a 3-D system, we 
use the well known relations [70] for the Fermi level energy Ep ~ h'^/Smd'^ = Eg and 
ek = ^NEp and obtain 



^fermi ^ ^ ^^^^ ^ 1.157V£„. (31) 



Howver, for 1-D system, we find ex ~ \NEp and Ep — to get 



Efermi ^ ^ ^NEp ^ NSo + ^NEo ^ IMNSo. (32) 

From Eqns. 31 and 32, it is clear that E^^^"^^ (both for 3-D system and 1-D system) is 
only marginally higher than E^°^^ (Eqn. 29) which indicates that X-motions have very 
little contribution to the G-state of a fermionic system (a; 15% in case of 3-D system, cf. 
Eqn. 31 and pa 1% in case of 1-D system, cf. Eqn. 32). 
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(v) . Process of Reaching G-State: In what follows from the above discussion the free 
energy (F) of the system can be identified to be a sum of two contributions F{K) and 

respectively, arising from K— and q— motions of particles. While F{K) refers to 
a gas of non-interacting quantum quasi-particles of mass 4m constrained to move in one 
direction only, F{q) could be equated, to a good approximation, to NEo as shown for 
the electron fluid ( a fermionic system ) in a conductor [66] as well as for liquid '^He ( a 
bosonic system ) [64]. As concluded by Mermin and Wagner [71], we note that F{K) is not 
expected to show any transition in case of a 1-D system, while F{q) — Nso does not depend 
explicitly on T and P. Consequently, there is no means to find a relation for a T and P 
for any possible transition in our system except by using possible physical arguments. In 
this context, we find that for a given particle density of the system, {d — A/2) decreases 
with decreasing T and at certain T — Tc, when [d — A/2) (with d — L/N) vanishes at 
large, q motions get frozen onto zero point motions oi q — qo — Tr/d and the system moves 
from a state of A/2 < o? to that of A/2 = d. While A/2 < d corresponds to 0(= 2qd) > 2n 
which represents randomness of (p positions and, therefore, a disordered state in 0-space, 
X/2 = d defines a state of = 2tt,- an ordered state in 0-space. Thus the system, on its 
coohng through Tc, moves from its disordered state to ordered state of its particles in the 
(p— space but this change does not represent an onset of a finite T phase transition. This 
is because the change refers to the relative configuration of particles (relative position 
(x = d), relative momentum [k = 2q) and relative phase position (0 = 2qd)) which 
remains unaltered with fall in T from T = Tc to T = 0. Evidently, the said change at 
T — Tc leads all the N particles to assume a configuration of the ground state of their 
0'— motions which naturally have no thermal energy at all T < Tc. In other words, the 
relative configuration of the system is effectively at T = and the thermal energy of the 
system at all T < Tc comes entirely from the i^— motions only. This clearly means that 
the said order-disorder of particles is effectively a T = change which does not violate 
Mermin- Wagner theorem [71] forbidding a transition in 1-D and 2-D systems at finite T. 

(vi) . Merger of N\ micro-states : We note that with all = go, different ^^^^ of 
(Eqn. 25) become identical and the latter attains the form of a single ^„ (Eqn. 22) 
which implies that all the A^! microstates merge into one at T = Tc and the entire system 
attains a kind of oneness, as envisaged by Toubes [72] ; the system at T < Tc is therefore 
described by 

*n(5) = (P^iq^mK), (33) 

where (pni^o) and (/)^(A'), respectively, define the relative configuration of the ground state 
and the collective excitations of the system. 

(vii) . Quantum Correlation Potential : The inter-particle quantum correlation potential 
(QCP) originating from the wave nature of the particles can be obtained by comparing the 
partition function (under quantum limit of the system), Zg = J2n^^p{~Fn/ kBT)\^niS)\'^ 
and its classical equivalent, Zc = J2n^w{~Fn/ ^bT) exp{—Un/kBT). Here ^n{S) is given 
by Eqn. 33. This can be justified because $„(S') is basically nothing but a product of 
paired wave functions obtained by simple superposition of plane waves and our conclu- 
sion that two particles in their SMW configuration satisfy d > \/2 screens out the HC 
interaction. Simplifying Un one easily finds that the pairwise QCP has two components. 
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A Ufj, pertaining to the g-motion of particles, controls the 0— position of a particle and 
we have 

^ = -ksTo ln[2 sm^ {kijXij/2)] , (34) 

with kij = kj — kj and Xij = x^—xj and T replaced by because T equivalent of q motion 
energy at all T < Tc is Tg. Ufj has its minimum value (—/csTo In 2) a,t (j) — {2n + l)7r 
and maximum value (oo) occurring periodically at A0 = 2n7r (with n — 1, 2, 3, ....). We 
note that Ufj increases by 

lc(50)2 = ^kBn{5<f>)' (35) 

for any small change (50) in around (f> — (2n + l)7r where it has its minimum value 
(—/c^To In 2). This indicates that particles experience a force = -CS(j) (force constant C 
= ^ksTo) which is naturally responsible for S(f) = and an ordered state in (/)-space. 
However, since Ufj is not the real interaction that may manipulate d, the state of order is 
achieved by driving all q towards qo- 

The second component of QCP, pertains to plane wave K motions, and it can be 
expressed by 

/ 27r|X2-Xi| 
1 ± exp 



Uij = —ksT In 



(36) 



by following standard procedure [73 and 74] applied to plane wave motion of non- inter- 



acting particles; while At = h/ ^2Ti{Am)kBT representing the thermal de Broglie wave 
length is associated with i^— motions (a kind of free motion where each particle appears 
to have 4m mass), +ve{—ve) sign stands for a bosonic(fcrmionic) system. In case of a 
bosonic system, Uij has its minimum (= —ksT \b.2) at \X2 — Xi\ = implying that Uij 
facilitates bosons to occupy a common X point. However, Uij for a fermionic system has 
its maximum (= oo) at X2 — Xi and, therefore, forbids two fermions from occupying a 
common X. 



(viii). Negative Thermal Expansion : In what follows from Sections-4(vii), the system on 
its cooling moves from its dis-ordered state in 0— space at T > Tc to an ordered state at 
T < Tc. The operational force for this change lies with : (i) QCP [U^j, Eqn. 34] which 
modulates the 0-positions of particles at A0 = 2n7r, and (ii) the zero point repulsion 
(c/. Section-3.(vi)) which keeps two particles at < x >> A/2. Consequently, when A/2 
increases with falling T, each particle pushes its neighbors to make space for its increased 
quantum size (A/2). However, the increase in A/2 virtually stops at Tc when g'-motions of 
particles get frozen at g = go (or A/2 = d) because the box gets totally occupied by closely 
packed wave packets of particles with L = NX/2. Evidently the zero-point repulsion of 
a particle on its nearest neighbor, at this stage, renders a force which tries to increase the 
box size, L. In all practical situations where forces restoring L are not infinitely strong, 
a non-zero increase {+SL) in L leading to a —ve thermal expansion coefficient of the 
system is expected around Tc. The experimental observation of —ve thermal expansion 
co-efficient would, obviously, provide proof for our predictions of : (i) the freezing of 
g-motions at go and (ii) the onset of an order in the 0— positions of particles. In this 
context, it may be noted that in agreement with our similar prediction for a 3-D system. 
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liquids "^He and ^He are really found to have -ve thermal expansion coefficient around 
2.17K and 0.55K, respectively [75]. 

(ix). Estimation ofTf. : We note that the potential energy contribution to net energy of 
the system is zero, since < A6{xij) > vanishes for all pairs of particles. Even zero-point 
energy So = h'^/SmcP, seemingly having potential energy character as indicated by its d 
dependence, represents the energy of residual g-motion which implies that the net energy 
of the system is kinetic. Consequently, to a good approximation, the lower bound of Tc 
could be equated to Tqse (the T equivalent of the G-state energy of a particle) with 
Tgse = El^'^/N (Eqn. 29) for the bosonic system and Tgse = E^^'^'^yN (Eqn. 30) for a 
fermionic system. However, to obtain more accurate Tc we need to account for T^^ = E^^, 
-the energy of thermal excitations (representing the correlated motions of particles) 
present in the system at Tqse and this renders 



T,^Tgse + T, 



(37) 



To estimate Tg^ for a bosonic system we note that the freezing of g-motions at g = g„, im- 
mediately follows an onset of the condensation of particles in = state because at this 
stage the system only has motions that may lose energy with faUing T. Although, as 
concluded in Section-3(v), the change at Tc is an effectively T = transition but one also 
finds that a quasi 1-D system realized in different laboratories basically represents a 3-D 
system where two dimensions are mechanically reduced to the order of inter-atomic sep- 
aration by increasing two corresponding components of momentum of confined particles 
and this agrees with the fact that these quasi 1-D systems do exhibit EEC in a manner 
identical to 3-D systems. Evidently, the K = condensation in our quasi 1-D laboratory 
systems can be identified as the BEG of non-interacting bosons because particles for their 
X-motions arc represented by plane waves. Since the particles for the /T-motions appear 
to have 4m mass, we can replace T^^ by \.Tbec to obtain 



1 /)2 

rpbose I 



1 

55 + 



(—) 



2/3' 



(38) 



for a bosonic system, where particles are free to move within its volume V on a surface 
of a constant potential. Similarly, for a system of bosons confined to a harmonic trap we 
have 



-ibose 



To + tT, 



BEC 



(39) 



where we use our result To = huoN^'^ / {nkB) [65] as well as Tbec = hw{N / 1.2^2)^'^ /kB 
obtained for such systems by Groot et. al. [76]. 

The Tj'°'"^ of a hypothetical 1-D system of HG bosons may be equated to To as its lower 
bound and to T^^"*"^ value given by Eqn. 38/39 as its upper bound since the /^-motion 
energy in such a 1-D system is expected to be lower than that present in a 3-D system. 
This indicates that Tg which represents about 66% of T^"^*^ in Eqn. 38 and about 60% in 
Eqn. 39 would be more close to T^°^'^ of a hypothetical 1-D system. 

In order to find T/'^''™* (the Tc for a Fermi system), we similarly use Eqn. 30 to find 
Tgse{= E^'^^^^^N) as its lower bound and assess Eex at Tgse to find equivalent Tex- 
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We note that such Eex should be a small fraction of E^^^"^^ because at such a low T 
not many fermions are in excited state and y/^''"** can be placed close to T^se ~ 
(a T slightly above To). It is important to note that T/^''"*' = T+ represents a point at 
which (i) g-motions get frozen at q = Qo, (ii) particles define an ordered state in the phase 
space and (iii) the system exhibits —ve thermal expansion coefficient. This should not 
be confused with a point below which a Fermi system transforms into a quantum fluid. 
As analyzed and discussed in [65, 66], the superfluid transition point in a Fermi system 
falls much below T^. This is because the fermions always have non-zero K and according 
to Pauli exclusion principle two fermions with equal q have to have unequal K or vice 
versa. Evidently, when they have equal they have unequal q and this possibility does 
not allow the relative configuration to stabilize with A0 = 2n7r unless the thermal energy 
falls below their binding energy which can be there only if the particles have inherent or 
induced inter-particle attraction [66]. 

5. Concluding Remarks 

(i) . The paper uses a new scheme to analyze the G-state of a 1-D system of N HC quantum 
particles. It concludes that each particle in the system should better be described by a 
macro-orbital (a kind of pair wave form) since each particle represents a part of the pair of 
particles {the basic unit of the system) having equal and opposite momenta (g, -q) at their 
CM which moves with momentum K in the laboratory frame. This agrees with the fact 
that particles in the system interact through a two body potential and their individual 
momenta (kj with i = 1,2,3, ...N) do not define good quantum numbers as soon as the 
particles have their wave mechanical superposition. Moreover the fact, that a particle in 
its macro-orbital representation has two motions [viz., the q— and K— motions] of the 
pair, -it represents, is consistent with Eqn. 6. 

(ii) . While all particles in the G-state of both systems (bosonic as well as fermionic) retain 
g-motions of identically equal q = qo = Tr/d, they all have K — for bosonic system and 
follow Fermi-distribution to occupy different states {K — nn/L with n — 1,2,3, ...) 
with K ranging from tt/L fa to Kp in case of a fermionic system. 

(iii) . Particles in the G-state of both systems define a close packed arrangement of their 
representative wave packets (each of size A/2 = d) with inter-particle phase separation 

= 2n7r (n = 1,2,3, ...), nearest neighbor distance < x >= d, and this packing does not 
allow them to have relative motion. Consequently, the G-state is a state of collisionless 
motion which can be identified as the residual zero-point g-motion. It differs from the 
higher energy states where particles can be perceived to have collisions because their 
effective size (A/2) is smaller than the average space {d) available to each particle. 

(iv) . While a bosonic system is expected to exhibit —ve thermal expansion coefficient at 
T xbose^ ^ fermiouic system should show such effect around T/^**™ pa Tq (c/. Section- 
4(viii)). This agrees with similar conclusion for a single particle in 1-D box [61] and two 
HC particles in similar box [58]. Our scheme applied to 3-D systems [63-66] also predicts 
similar expansion whose accuracy is established by experimentally observed expansion of 
systems like liquids '^He and ^He [75]. 
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(v) . The G-state configuration of our system (c/. point (ii), above) does not assume 
stability against any perturbation {viz., the flow of the system) because the particles 
in this configuration are in a state of persisting zero-point repulsion and the necessary 
inter-particle attraction to counter this repulsion is absent. If such an attraction exists 
and the system remains fluid even at T < T^, then the system can, in principle, exhibit 
phenomenon like superfluidity below certain T (say Tx); our analysis of 3-D systems 
such as hquid ^He [64] and liquid ^He [65, 66] concludes that Tx for a bosonic system 
falls close to To but for a fermionic system it is expected to be orders of magnitude 
lower than To because X— motion energy retained by fermions due to Pauli exclusion is 
good enough to de-stabilise the G-state configuration unless kBT falls below per particle 
binding energy. In our forthcoming paper, we plan to study different aspects of our 1- 
D system when weak interparticle attraction is added as a perturbation. However, our 
detailed study of a 3-D bosonic system [64] and a brief quahtative analysis of widely 
different many body systems (including low dimensional systems and fermionic systems) 
[65] clearly reveals that a 1-D system realised in a laboratory must have basic properties 
such as superfluidity and related behavior of 3-D systems if weak inter-particle attraction 
is present. As shown in Section-4(v) the predicted order-disorder transformation of our 
1-D systems is consistent with Mermin Wagner theorem [71]. In this context, it may be 
emphasized that the transformation is basically a quantum transition expected to occur 
at T = 0; however, it occurs at a non-zero T for the proximity of zero-point g— motions 
{q = Qo) with i^— motions which solely represent the thermal motions of particles at 
T<To. 

(vi) . The formation of a SMW from the superposition of two plane waves of two HC 
particles is as natural as the phenomena of interference and diffraction of particles such 
as strongly interacting electrons, neutrons, helium atoms, etc. [77, 78]. Since the nature of 
interference and diffraction patterns of these strongly interacting particles does not differ 
from the nature of such patterns for non- interacting photons, it is evident that only 
wave nature (not the inter-particle interactions) modulates the relative phase positions of 
particles in their wave superposition. Evidently, the fact that these experiments support 
the formation of a SMW can not be doubted. 

(vii) . We understand that, in principle, two particles described by plane waves should 
have their superposition independent of their d and A. However the experimental fact, 
that wave nature dominates the behavior of a many body system only when X> d [59, 69], 
indicates that effective wave superposition of two particles becomes possible only when 
X> d and this could be so because no particle in nature manifests a real plane wave. Since 
the important consequences of wave superposition (c/. Section-3.0) such as zero point 
repulsion leading to < x >> A/2, inter-particle 0— correlation which arranges particles 
at A(j) = 2mr and renders coherence of particle motion, etc. (incorporated in the present 
analysis) are not used in [1-6] , our results are expected to differ from those of [1-6] . Guided 
by this observation and the experimental support to the formation of a SMW (Point-vi, 
above), we may highlight the said differences by comparing our results with those of [4] for 
7 = cxD case. If our inference for allowed q = qo = ir/d and = tt/L ^ for the G-state of 
a bosonic system is used to determine the corresponding momenta through kj = q + K/2 
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and kj = —q + K/2 of two particles in superposition, each boson should either have 
7r/d or —n/d momentum and we find that these values differ from possible k = ±sTr/L 
(with s — 1,2,3, ...) considered in [4]. Similar estimates of k for a fermionic system using 
q = Qo = Tr/d and K ranging between n/L ^ and Kp = Nn/2L = n/2d render a band of 
allowed k ranging from —q + Kp/2 = Sir/Ad to q + Kp/2 = bn/Ad (with least difference of 
ir/iL) which do not match with k = ±2stt/L considered in [6]. Note that these allowed k 
have been determined just for their comparison with those of [1-6] , otherwise our analysis 
reveals that k, and k^ do not represent good quantum numbers for the two particles in a 
state of wave superposition. Similarly, per particle G-state energy So = h'^/Smd,'^ obtained 
by us for a bosonic system is three times higher than (l/3)£o conchidcd in [4]. It may 
be noted that our results reveal a well defined picture of the G-state of 1-D systems. 
Accordingly, the system in its G-state is a close packed arrangement of the wave packets 
of its particles with < x >— d, q — qo and = 2mr and this fact has been revealed 
for the first time. One may find a number of other points where our results differ from 
those of [1-6]. However, since an effective wave superposition of particles docs not exist 
in higher energy states of X < d, such particles can equally well be described by plane 
waves as considered in [1-6] . Evidently, such studies seem to provide reasonably accurate 
understanding of higher energy states (or high T phase) of the system and the present 
study, revealing the low energy states (or low T phase), can be identified to supplement 
the results of [1-6] in providing a complete understanding of a 1-D system. In view of an 
equivalence in the dynamical behavior of two HC particles interacting through Vhc and 
AS{x), particularly, for low energy states of A/2 > cr, it can be stated that our analysis 
also applies to particles of finite a satisfying A/2 > a. Finally, as an important inference 
of this study, the true picture of the low energy states of our system can be revealed if 
its theory incorporates the impacts of the wave superposition and zero-point repulsion 
{obvious effects of wave particle duality), viz., < x >> A/2 q > n/d and A0 > 2n7r on 
the relative configuration of two particles. 
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